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Abstract. Near a black hole or an ultracompact star, motion of particles is 
governed by strong gravitational field. Electrically charged particles feel also 
electromagnetic force arising due to currents inside the star or plasma circling 
around. We study a possibility that the interplay between gravitational and 
electromagnetic action may allow for stable, energetically bound off-equatorial 
motion of charged particles. This would represent well-known generalized 
Stormer's 'halo orbits', which have been discussed in connection with the motion 
of dust grains in planetary magnctospheres. We demonstrate that such orbits exist 
and can be astrophysically relevant when a compact star or a black hole is endowed 
with a dipole-type magnetic field. In the case of Kerr-Newman solution, numerical 
analysis shows that the mutually connected gravitational and electromagnetic 
fields do not allow existence of stable halo orbits above the outer horizon of black 
holes. Such orbits are either hidden under the inner black-hole horizon, or they 
require the presence of a naked singularity. 
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1. Introduction 

Investigation of charged particle motion in strong gravitational and electromagnetic 
fields ranks among the elementary exercises in theoretical physics and astrophysics 
related to black holes or compact stars [1-6]. It helps us to reveal the structure of 
force fields acting on the particles. Numerical integration of the motion equations is 
very helpful, but often fails to give full understanding. On the other hand, analytical 
techniques are often restricted, for practical reasons, to the motion along symmetry 
axis [7], and to the equatorial [8] or spherical motion [9], or to the special case of 
free-infalling particles [10]. The general off-equatorial motion in the field of rotating 
black holes or compact stars still remains an interesting open problem deserving 
attention. In this paper, we tackle this problem focusing on the case of halo orbits, 
i.e., off-equatorial energetically bound orbits that are stable with respect to small 
perturbations, but they never cross equatorial plane. 

As a motivation of our work, we mention one of the milestones of the classical 
space physics in this branch, i.e., Stormer's analysis of charged particles motion in 
a purely dipolc magnetic field (classical Stormer problem) [11], providing us with 
the basic physical description of radiation belts surrounding a magnetized planet. 
Radiation belts are known to be composed of individual ions and electrons whose 
motion is governed by magnetic forces only. In order to describe the dynamics of 
charged dust grains in planetary magnetospheres, when there are much smaller charge 
to mass ratios, one can employ the 'generalized Stormer problem', where both the 
planetary gravity and co-rotational electric field play a role [12]. Such studies point 
out the existence of the dust grains halo orbits near, e.g., Saturn [13]. 

Could such halo orbits survive in strong gravitational fields near compact objects 
with additional magnetic or electric fields? The situation is not obvious in case of exact 
solutions of Einstein-Maxwell equations, where the gravitational and electromagnetic 
fields are mutually interconnected. Surprisingly enough, this problem has not yet been 
fully answered in the case of well-known Kerr-Newman black hole. In fact, we can 
anticipate such orbits to occur near the black hole horizon, where gravity is strong. 
Their exact location should depend on the type of particles we consider, namely, on the 
charge to mass ratio. Halo orbits could have outstanding astrophysical consequences 
for radiation of accreting compact objects. This is because the halo orbits form two 
lobes of stable motion - one above and the other one below equatorial plane - through 
which the accretion disc radiation has to pass on its way towards a distant observer. 
Therefore, the halo lobes act as a part of radiation reprocessing corona. 

We investigate the problem with three qualitatively different kinds of compact 
objects. These are magnetized compact star, studied in terms of pseudo-Newtonian 
model (section |2|, Schwarzschild black hole with a current loop (section [3]), and Kerr- 
Newman black hole or naked singularity (section The first case is most close to 
the original Stormer's approach, except that we consider strong gravity. The second 
case may serve as a toy model of the spacetime associated with a magnetized disc 
about a Schwarzschild black hole. In the last case, the solution was partly found 
in [14,15], where the authors conclude that the halo orbits do exist in Kerr-Newman 
spacetimes. Nevertheless, because of many input parameters characterizing the central 
object and the particle motion itself, the detailed discussion of the orbit existence was 
not presented there. We can remark that the rotating black hole immersed in a 
magnetic field will establish a stationary system with a non-zero net charge [16]. In 
astrophysically realistic magnetic fields, this charge is probably very small. However, 
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we adopted general approach to the problem, which allows us to explore the whole 
parameter space in a systematic way. Thus, we considered situations in which the 
charge is small, but also those with large charge for completeness of our discussion. 

The paper [17], dealing with a similar problem, i.e., with the stability of charged 
spherical shells in Kerr-Newman spacetimes, suggests that the existence of stable halo 
orbits is restricted to the naked-singularity spacetimes only. The authors of that paper 
employed the traditional four-dimensional covariant form of the motion equations. 
Our complete discussion of the Kerr-Newman spacetimes is based on inertial forces 
analysis [18-21], which, in combination with the effective potential, gives a semi- 
analytic condition for the halo orbits existence. 



2. Compact star with dipole magnetic field 

We consider a rotating object of mass M, radius R and spin f2 oriented in positive 
^-direction (z+). The object is a source of central gravitational field (potential tp) 
described in terms of Gaussian units and Cartesian coordinates (x, y, z). Further the 
object is endowed with rigidly co-rotating aligned dipole magnetic field B — V x 
Magnitude of the magnetic field measured at the object equator is B , and the induced 
co-rotating electric field has intensity E — ~ ^(Sl x r) x B, where r = (x, y, z). 

Similar set up of the problem was employed in [12], but here, we describe 
gravitational field of the central object by the Paczynski-Wiita potential [22] 

/ - GM m 
r - r s 

where rg = 2GM/c 2 is the Schwarzschild radius. Further, we use the common form 
of the dipole vector potential 

M 

A=-f(-y,x,0), (2) 

where M. = BqR 3 is magnitude of the z+ oriented dipole magnetic moment. 
The particle motion is described by Hamiltonian 

where 11 is the potential energy, and p, m, q are, respectively, the momentum, mass and 
charge of the particle. In the cylindrical coordinates (p, <j>, z), where p — (x 2 + y 2 ) 1 ^ 2 
and <p — arctan(2/ / x) , 

n = ^(p 2 p + P 2 z ) + u cS , (4) 

with the effective potential 

U eS = ^ m -^(p ( j > --y^) 2 + a s milj + a^m. (5) 

In the above equation, p^ — mp 2 d<f>/dt + 7^ is the conserved angular momentum, 
7 = qAi/c, and ^ = p 2 /r 3 is the dipole stream function. This function enables us to 
express the electric field in the form 

E = --MttW. (6) 

c 



X Electric and magnetic fields are related by the force-free relation E + (l/c)V X B = 0, where 
V = n x r. 
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Figure 1. Effective potential U a g and its contours for motion of particle with 
(in Gaussian units) (a) m = 1.67 X 10" 24 , q = 4.85 X lCT 10 , = -3 X 10~ 8 
(counter-rotating proton), or equivalently q = —4.85 X 10 — 10 , = 3 X 10 — 8 (co- 
rotating hydrogen anion); (b) m = 4.2 X 10 -12 , q = —3.3 X 10~ 5 , = 2.5 X 10 4 
(co-rotating negatively charged dust grain), in central gravitational and rotating 
dipole magnetic fields, corresponding in the zero approximation to the fields 
of magnetized compact star with M = l.SAfg, R = 10 6 cm, and (a) £7 = 0, 
B = 10 3 G; (b) H = 10 3 rad/s, B = 10 8 G. Positions of stable halo orbits 
correspond to the potential minima. Corresponding, stable in radial direction, 
circular equatorial orbit corresponds to the saddle point in the equatorial plane. 

Following [12], we also introduce 'field switches' <r g and <j r , taking values or 1. 
Characteristic behaviour of the effective potential ^ is shown in figures [l] and [2] and 
suggests the existence of halo orbits near magnetized compact stars. 

Measuring time by 2-7r/f2, distance by 1Z — (GM/fi 2 ) 1 / 3 , and scaling the mass by 
the particle mass m, we find 

H =l(p 2 p +f z ) + U eS , (7) 

U eS = Ul- 5 lL) _ ff l_ +(Tr ^! ( 8 ) 
2 \P r J r — rs f 

where the hat denotes scaled quantities, and 

. _ ns . 2GM n 7 

Introducing the scaled spherical coordinates (r,<f>,9), where f — (p 2 + z 2 ) 1 / 2 and 
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Figure 2. Effective potential U a g and its contours for motion of particle with 
(in Gaussian units) (a) m = 4.2 X 10~ 12 , q = 3.3 X 10~ 5 , = — 5 X 10 4 (counter- 
rotating positively charged dust grain) (b) m = 4.2 X 10~ 12 , q = 3.3 X 10~ 5 , 
f>0 = 10 X 10 4 (co-rotating positively charged dust grain), in central gravitational 
and rotating dipolc magnetic fields, corresponding in the zero approximation 
to the fields of magnetized compact star with M = 1.5Mg, R = 10 6 cm and 
f2 = 7 X 10 3 rad/s, Bq = 10 8 G. Positions of stable halo orbits correspond to the 
potential minima. Concomitant unstable circular orbit corresponds to the local 
maximum in the case (b). 



9 = arccos(z/f), and the orbital angular velocity of a particle 

dt r 2 sin 9 r s 

we can write the effective potential in the form 

1-2-2-2/! °~k a Y 8 sin 2 9 

U eS = -u?r 2 sin 2 9 - -^A- + ; • 11 

2 r — rg r 

Conditions for the existence of the halo orbits, dfU c s = and deU c s — 0, are: 

-9- . 9o 5(cb — o- r ) sin 2 9 a„ „ 

-Ld 2 f sin 2 9+ -i £ + g ., ? =0, 12 

r z (r — rs) z 

-I cos 9 sin <9(w 2 f 3 + 2uS - 2cr r (5) = 0. (13) 
r 

Eliminating Co from these equations, we can write the halo orbits existence condition 
in the form 

±a 2 f 3 esc 6 

5= I . (14) 

cr g (f — rs)[(6(T g f) 1 / 2 + 3<7 r r(r — fs) sin 6*] 
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Thus, for arbitrarily fixed coordinates of halo orbits, there are real values of the 
parameter 5, which can be set by tuning the object and particle parameters. From 



equation ( 14 1, it is clear that the co-rotational electric field is not necessary for 
the existence of the halo orbits, i.e., the rotational effect of the magnetic field is 
not essential, although it influences the properties of the halo orbits. On the other 
hand, there are no halo orbits without the gravitational field. These conclusions are 
qualitatively the same as in [12]. 

In the presence of the gravitational field and absence of the electric field, there are 
co-rotating stable halo orbits for particles with negative charge and counter-rotating 
orbits for particles with positive charge. The effective potential is invariant under 
simultaneous charge reversal and angular momentum reversal. In the case of presence 
of both the electric and magnetic fields, i.e., adding the rotation to the magnetic field, 
there are no qualitative changes for negatively charged particles motion along (always 
stable) halo orbits. Such particles must co-rotate with the source. On the other hand, 
the positively charged particles can either counter-rotate with the source or co-rotate, 
moving along stable halo orbits (§] 

3. Schwarzschild black hole with dipole magnetic field 

We consider a Schwarzschild black hole of mass M described by the metric [24] 

ds 2 = - (i - ^f) dt2 + (l - ^) 1 dr2 + r2 ( d ° 2 + sin2 ed ^XV>) 

in the geometric units (c = G = 1) and standard Schwarzschild coordinates (t, r, <p, 6). 
The black hole is surrounded, in the equatorial plane, by a toroidal current loop of 
radius R, which is supposed to contribute to the electromagnetic field but not to the 
gravitational field. The magnetic field generated by the loop is influenced by the black 
hole gravitation [25]. It can be properly described by the leading (dipole) term of the 
full multipole solution [26], i.e., by 



Ai = --6. 



3 ^Mr 2 sin 2 6 



^ ! 2M\ . 2M . 2M 2 1 



n = -g ij (P i -qA i )(p j -qA j ) > (17) 



M3 > > > < 16 > 

in the region r > R. Here Ai — nR 2 I(l — 2M/r) 1 / 2 is magnitude of the dipole moment 
of the magnetic field and / is the current in the loop. 

The 'super Hamiltonian' for the motion of a test particle reads 
1 

2 J ^ 

where pi is the generalized momentum. The first Hamilton's equation implies 

dx l /dX — p 1 — qA l = p 1 . The second Hamilton's equation ensures that the momenta 

/ 2M\ dt , , 

p t =pt + qA t = -m\l-— J — = -E (18) 

P<j>=P<t> + lAj, = mr 2 sin 2 9— + qA^ = L (19) 

are conserved. Introducing now the specific energy, angular momentum and charge 
by the relations E = E/m, L = L/m and q = q/m, we find 

E 2 = (f) 2 + r 2 (l- 2 -^)(p 9 ) 2 + V 2 s . (20) 



§ Preliminary analysis of differences between the original Newtonian (weak gravity of a planet) and 
pseudo-Newtonian (strong gravity of a compact star) was presented in [23]. 
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Figure 3. Effective potential V c g and its contours for motion of charged particle 
with L = 1 (in geometric units and units of M) in the Schwarzschild spacetime 
with the test dipole magnetic field of current loop characterized by the parameter 
k = —4. Positions of stable halo orbits correspond to the potential minima. 
Concomitant equatorial circular orbit, stable in radial direction, corresponds to 
the saddle point in the equatorial plane. The effective potential is not relevant 
in the region limited from above by the radius of the current loop R = 6 (light 
gray), entirely containing the dynamic region r < 2 of the spacetime (dark gray). 



The quantity V e g is the effective potential taking in the dimensionless coordinates 
r/M — > r, t/M — > t the form 

2 ' 



Kff = (!--)< 1 + 



L 



r sin f 



kr sin 8 



In 1 



where 



3, „ 2r / 2M\ 1/2 

k = -qnR III = const. 

8 \ r J 



(22) 



Figure [3] shows the characteristic behaviour of the effective potential (21), suggesting 
the existence of halo orbits. 

In order to confirm the presence of halo orbits, one needs to check conditions 
for the existence of halo orbits d r V c fi = and deV c ff = 0, which now imply the 
'master' relation for the parameter k 



±ircsc( 



-[2 + 2r + r 2 ln(l - ^)][6 + 2r + r 2 In (1 



(23) 



The reality condition in this relation restricts the halo orbits existence to the region 
r > 2.3231. But this condition is automatically satisfied due to the assumption of the 
loop radius R > 6 [26]. Thus, considering a circular orbit with the radius r > R and 
latitude < 9 < 7T, we can find corresponding value of the parameter k and fit this 
value by tuning parameters of the particle, black-hole and magnetic field in equation 
(22 1. The Schwarzschild black hole with a test dipole magnetic field and a non- 



rotating compact star with aligned (static) dipole magnetic field can be considered as 
equivalent to each other. As intuitively expected, there is the same type of behaviour 
of the relativistic effective potential ( [21) and the pseudo-Newtonian one §5§ without 
the contribution of the electric field. This consistency in type of the effective potentials 
behaviour confirms our results from section[2] The other types of the effective potential 
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(|5| behaviour are not recovered in the Schwarzschild case, simply because there are 
no rotational effects inducing the electric counterpart. 

4. Kerr-Newman spacetime 

Now, we study halo orbits existence within Kerr-Newman solution of the full Einstein- 
Maxwell equations. In the dimensionless Boyer-Linquist coordinates (t, r, </>, 9), the line 
element takes the form [24] 

A S 
ds 2 = - —(dt - asm6d<j)) 2 + — dr 2 + Ed9 2 

+ S ^-[(r 2 + a 2 )d^-adt] 2 , (24) 

where A = r 2 — 2r+a 2 +e 2 and £ = r 2 +a 2 sin 2 6. Quantities a and e are the rotational 
and charge parameters of the spacetimes. Non-zero components of the antisymmetric 
electromagnetic field tensor F^j = Ajj — Aij read 

F rt = e(r2 -g C ° s2g) , (25) 
— ae sin 2 9{r 2 — a 2 cos 2 9) 

trit) = ^ ' ( 26 ) 

-q 2 ersin26> 
*et ^ ' 

aer sin 26>(r 2 + a 2 ) 
= ^ • ( 28 J 

The effective potential for the motion of a particle with the specific charge 
q = q/m, takes the form [24,27] 

W«= Y -±^EM, (29) 

where 

X= (r 2 + a 2 ) 2 -Aa 2 sm 2 9, (30) 

y = (La + qer)(r 2 + a 2 ) — LaA, (31) 

Z = (La + qer) 2 — AE — AL 2 / sin 2 6», (32) 

and L — L/m is the conserved axial component of specific angular momentum. 



Characteristic behaviour of the effective potential ( 29 1 suggesting the existence of 
halo and equatorial orbits in Kerr-Newman spacetimes is shown in figures |4j|6] 

In order to study the halo orbits, we find it very helpful to employ the 'force 
approach' [28] [jj] For uniform circular motion of a particle at constant radius and 
latitude, and with the velocity v measured by ZAMO [29], we can write (see Appendix 
A) two force equations: 

- Qr - (lv) 2 Z r — j 2 vC r — q^(£ r + vM r ), (33) 

-Go - (ivfZe - i 2 vC 9 = qi(£ 9 + vM e ), (34) 

where Lorentz factor 7 = (1 — v 2 )" 1 ! 2 . Here Qj, Zj, Cj denote the mass and velocity 
independent parts of the gravitational, centrifugal and Coriolis inertial forces, and 

|| We found [21] the force analysis of circular geodesies to be much more efficient and straightforward 
in comparison with the effective potential approach. 
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Figure 4. Effective potential W c s and its contours for motion of charged 
particles with L = 1 and q = — 5 (in geometric units and units of M) in the 
Kerr-Newman (a) black-hole spacetime with a 2 = 0.04 and e 2 = 0.95; (b) naked- 
singularity spacetime with a 2 = 0.6 and e 2 = 0.95. Positions of stable halo orbits 
correspond to the potential minima. Concomitant equatorial circular orbit, stable 
in radial direction, corresponds to the saddle point in the equatorial plane. The 
effective potential W B ff is not relevant in the region between the event horizons 
(gray), where A < 0. 



£ j and Mj stand for the charge and velocity independent parts of the electric and 
magnetic forces (see relations ( A. 19 1— ( A. 30 1). 



Eliminating q from equations (33 1 and ( 34 1 , and assuming < 8 < n/2, i.e., 
omitting stationary equilibrium points on the axis of symmetry and circular orbits in 
the equatorial plane, we obtain a cubic equation 

Av 3 + Bv 2 



Cv + D = 0, 



(35) 



where 



A = Mg(G r - Z r ) + M r {Zg - g g ), 
B = £ r (G r - Zr) + £r(Z e - Go) + C e M r 
C = Cg£ r - C r £ 9 + g e M r - GrMe, 
D = £ r Ge — £eG r - 



C r M 6 



(36) 
(37) 
(38) 
(39) 

This gives three, in general complex, solutions vi(r,6;a,e), vjj(r,6; a,e) and 
vni(r, 6; a, e) for possible orbital velocities of charged particles moving along the halo 
orbits. Parameterizing them by the spin a and latitude 9, we can investigate the 




Figure 5. Effective potential W c f{ and its contours for motion of charged 
particles with L = 0.1 and q = 2.9 (in geometric units and units of M) in the Kerr- 
Newman black-hole spacetime with a 2 = 0.5 and e 2 = 0.3, in the regions (a) above 
the outer event horizon; (b) under the inner event horizon. Positions of unstable 
stable halo orbits correspond to the saddle points of the potential. Concomitant 
equatorial unstable circular orbit corresponds to the potential maximum in the 
equatorial plane. The effective potential W e ff is not relevant in the region between 
the event horizons (gray), where A < 0. 



validity of the condition v t 6 R and — 1 < < 1 (i = I, II, III), in the plane (r x e). 
Our numerical analysis of Uj confirms the existence of halo orbits in the Kerr-Newman 
naked-singularity spacetime as well as in the black-hole spacetimes under the inner 
event horizon and even over the outer horizon (see the representative case in figure [7]). 

We conclude our study of halo orbits existence by the discussion of their stability, 
namely by seeking for the stable ones. In the previous two cases (sections [2] and [3| , 
the stability of the halo orbits could be directly inferred from the behaviour of related 
effective potentials. For the present purpose, we can use equation (33 1 and express 



the specific charge qh of the particle moving along the expected halo orbit at velocity 
v by the relation 

g r {v 2 - 1) - v{C r + vZ r ) 



2h 



(£ r + vM r )Vl~ 



(40) 



Condition (40 1 provides three possible values qh,i related to Vi. By rewriting relation 



(19 1 for the Kerr-Newman case and using the 4- velocity decomposition (A.l), we can 
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Figure 6. Effective potential W e ff and its contours for motion of charged 
particles with L = 8 and q = —5 (in geometric units and units of M) in the 
region above the outer event horizon of the Kerr-Newman black-hole spacetime 
with a 2 = 0.2 and e 2 = 0.5. The potential minimum corresponds to the stable 
equatorial circular orbit. The effective potential W c h is not relevant in the region 
between the event horizons (gray), where A < 0. In all numerically tested cases, 
no stable halo orbits have been found, only the stable equatorial circular orbits 
were confirmed. 




(b) 



0.5 1.0 
ArcTan[r] 




0.5 1.0 
ArcTan[r] 



Figure 7. Existence of halo orbits for 6 = 1.1 in Kerr-Newman black-hole (BH) 
and naked-singularity (NS) spacetimes with a 2 = 0.4. Positions of horizons are 
denoted by the thick curve, (a) White areas correspond to the regions where the 
root vj is real and — 1 < vj < 1, i.e., where the halo orbits occur, (b) White areas 
correspond to the regions where the stationary points of the effective potential 
W e ff, corresponding to the values of vj, are minima, i.e., where the halo orbits 
are stable. 



write for the specific angular momentum of particles at halo orbits 

Lh = 7 u \/500 + 9hA/>, (41) 

obtaining three possible values L^.i related to gh,i and fi- 

The stable halo orbits must satisfy the conditions of effective potential minima, 

i.e., 

8 2 r W cS (r, 8; a,e,L = L h , q = q h ) > 0, (42) 
d 2 e W eS {r, 8; a,e,L = L h , q = q h ) > 0, (43) 
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which we analyze in a numerical way in the plane (r x e), considering the 
parameterization by a and 6, for all three pairs of L\ wi and q^. Searching 
systematically through the parameter space of (r x e x a x 6?), we found no combinations 
of parameters allowing the existence of stable halo orbits in the Kerr-Newman black- 
hole spacetimes above the outer event horizon. Representative results of such analysis 
are shown in figure [7J where both black holes and naked singularities are explored. 
In these plots we keep the spin a fixed, whereas the electric charge e varies along the 
ordinate. We find this type of plots very instructive, as they capture the entire range 
of radii from the origin up to spatial infinity. Shading distinguishes the areas where 
the halo orbits cannot exist (left panel), or where they are unstable (right panel). On 
the other hand, the region under the inner event horizon and the naked-singularity 
spacetimes do exhibit stable halo orbits. 

5. Discussion 

Our analysis has shown that the stable halo orbits can exist near rotating as well 
as non-rotating magnetized compact stars. Rotation enriches the classification of 
the stable halo orbits types. Our first example, the pseudo-Newtonian gravitational 
potential and the classical dipole magnetic field do not describe the situation with the 
same precision as a relevant solution of Einstein-Maxwell equations, especially in the 
region under the photon circular orbit, where the Paczynski-Wiita potential does not 
describe the gravitational field properly. However, the pseudo-Newtonian approach 
is sufficient to demonstrate how the halo orbits arise in strong gravity. The halo 
orbits can exist even far away from the photon circular orbit, for realistic values of the 
particle specific charge and the compact star magnetic field strength. Moreover, at 
such distances, the gravitational and magnetic fields of one of the investigated cases, 
i.e., the slowly rotating magnetized compact star, are nearly of the same character as 
the fields generated by Schwarzschild black hole with a current loop in the equatorial 
plane, our second example. The general relativistic study of such configuration has 
confirmed the existence of stable halo orbits. 

In the fields of Kerr-Newman black holes, the third example, no stable halo 
orbits have ever been, to our knowledge, found outside the outer horizon, although 
there is a magnetic field of dipole character as well. But here, in contrast to the 
case of magnetized compact stars and Schwarzschild black holes with the externally 
generated magnetic field, where the magnetic field strength is generally independent of 
the spacetime parameters, the magnetic (and related electric) field is determined by the 
internal spacetime parameters (the spin and electric charge) and cannot be imposed 
as an independent variable. It seems that the connection of the electromagnetic field 
with the spin of the spacetime prevents the existence of stable off-equatorial circular 
orbits outside the outer horizon. Only unstable halo orbits are allowed there. On the 
other hand, the combinations of spacetime parameters allow the existence of stable 
halo orbits under the inner horizon of black holes. The inner horizon is also a Cauchy 
horizon that may focus to a spacelike one [30] . This may imply that the continuation 
of the spacetime manifold beyond it lacks any physical or astronomical relevance. 
However, beyond the ring singularity of Kerr metric, there is an asymptotically flat 
universe. If we actually live in such an asymptotically flat region, in addition to other 
interesting phenomena like closed time like curves, one may also observe the halo orbits. 
Therefore, the existence and properties of halo orbits may be a useful tool probing the 
internal properties of black holes along with the traditional analysis of their geodesic 
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structure [31]. In principle, as follows from our analysis, halo orbits become accessible 
in the naked-singularity spacctimcs as well. Therefore this fact may be important in 
the identification of naked singularities. 

These conclusions are consistent with the results of a purely analytical study of the 
situation on the axis of symmetry, where the circular motion degenerates to a point. 
There, only unstable positions are allowed above the outer event horizon. The stable 
positions occur only under the inner horizon or in the naked-singularity spacetimes [7] . 
Thus, within the framework of the Kerr-Newman black-hole spacetimes, the only 
astrophysically relevant stable circular orbits at constant radius and latitude are those 
in the equatorial plane [32]. 

The equatorial plane is where a gaseous accretion disk can reside. Although 
the accretion discs are often described as geometrically thin and planar structures, 
in reality they must have extended atmospheres which reach far above and below the 
equatorial plane. This motivated us to speculate that the halo orbits could be relevant 
for accretion in strong gravity. High frequency kilohertz quasiperiodic oscillations 
(QPOs) observed in some microquasars and binary systems with compact stars are 
frequently explained by variety of models based on the equatorial quasicircular motion, 
with characteristic orbital Kcplcrian and epicyclic frequencies. Most promising seem 
to be the relativistic precession model [33], and the orbital resonance model [34,35] or 
its generalization to the orbital multiresonant model [36,37] for both binary systems 
with black holes (microquasars) and compact stars; in the case of near-extreme Kerr 
black-hole candidates (e.g. well known GRS 1915+105 microquasar) the complex 
high frequency QPOs patterns could be explained using the extended resonance 
model with the so-called hump-induced oscillations, additional to the orbital epicyclic 
oscillations [38,39]. By using the orbital resonance models, black hole parameters, 
especially the spin, could be determined [40-42]. 

The existence of two mutually detached off-equatorial lobes of bound stable orbits 
can have profound consequences for plasma oscillations near compact objects. It allows 
us to imagine a situation when the bulk motion is along the stable circular orbits in 
the center of the lobes, around which a fraction of trapped particles oscillate. We 
imagine that the oscillations modulate the electromagnetic signal from plasma cloud, 
and this could be observed as a mode of oscillations in the detected radiation. 

The halo orbits in the magnetic field of compact stars and black holes could be 
related to the oscillatory motion with 'halo' radial and vertical frequencies, which 
could be considered as complementary model of the QPOs observed in the binary 
system of compact stars and microquasars. 

6. Conclusions 

Our study has shown that a strong gravitational field endowed with a magnetic field 
of dipole character enable the existence of the so-called stable halo orbits, i.e., stable 
circular off-equatorial orbits. Such a composition of fields is expected in the vicinity 
of compact objects, e.g., in physically relevant situations with magnetized compact 
(neutron, quark) stars, black holes with current loops in the equatorial plane, or in the 
case of charged and rotating black holes. We have found that the stable halo orbits 
can take place near all of the three mentioned kinds of compact objects. But in the 
Kerr-Newman charged and rotating black- hole spacetimes, the stable halo orbits are at 
most of marginal astrophysical importance, being hidden under the event horizon. We 
mentioned the possibility that the observer could actually live in an asymptotically flat 
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region beyond the ring singularity, however, this would be clearly a highly speculative 
proposal. On the other hand, in the field of magnetized neutron stars or near black 
holes with current loops in the equatorial plane, the stable halo orbits appear outside 
the body and can be astrophysically relevant. 
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Appendix A. Projection of Lorentz equation and inertial forces formalism 

We start with the decomposition of the particle 4-velocity into the form [18] 

u l = -f(n l + vt 1 ), (A.l) 

where the 4-velocity field n l satisfies the conditions 

n k n k = -1, n l \7 t n k = V fc $, n [i V j n k] = 0, (A.2) 

and the vector t % is the unit spacelike vector orthogonal to it, along which the spatial 
3-velocity with magnitude v is aligned. The 4-velocity field n l can be chosen in the 
form 



„2* 



(A.4) 



-thi, (A.3) 

thus it corresponds to 4-velocity of stationary observers, parallel to a timelike vector 
field l\ 

Now, we can express the left-hand side of the Lorentz equation 
mu k V k Uj = qFjkU k 

in the form 

ma k = m[7 2 V fe $ + j 2 v(n l V z T k + t 1 V \n k ) 
+ 7 2 vVViT fc + (f7)r fe + jn k ], 
where (iry) = u l V ' i(^/v). By using the projection tensor 
h tk = g ik + nin k , 



we project the uniquely decomposed 4-force (A. 5) 

i , k 
nia^ = mhj a k 



-Gj 



z j Gj 



L 



j i 



(A.5) 
(A.6) 
(A.7) 



where 



Gj 
Z, 



m(jv) f*V»fj, 



(A.8) 
(A.9) 
(A.10) 
(All) 

can be interpreted as the gravitational, centrifugal, Coriolis and Euler inertial forces 
[18]. In the above relations Xj — n l (WiTj — VjTj) and V = — u l V i(u k u k v) . The vector 



Cj = — m7 V X 
Lj = - mVTj 



j ■ 
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f l = e*r l , with its covariant form fj = e~*Ti, is the spacelike unit vector parallel to 
r' in the so-called optical reference geometry, defined by the relation 

hik = e- 2 *h lk . (A.12) 



Projecting the right-hand side of the Lorentz equation (A.4), we obtain the 
equation [19] 

qh)F lk u k = E J +M J , (A. 13) 

where 

Ej =q<YF jk n j , (A. 14) 

Mj = qjv(F jk r k + n 3 F kl n k r l ) (A.15) 

are the electric and magnetic forces, respectively. 

We focus our attention to the case of the symmetric and stationary spacetimes 
with the Killing vector fields rf = SI and = S\, and consider n % to be 4- velocity of 
the ZAMO observers, i.e., 

n 1 =e"*(rf + ft Z AMof), (A.16) 

e 2 * = - (v l + OzamoC* ){m + ^ZAMoCi), (A. 17) 

where ^zamo = — ^'^i/C^Cr Considering uniform circular motion of particle at 
constant radius and constant latitude, i.e., 

r l = (£*&)" 1/2 £*. (A.18) 

in Kerr-Newman spacetime, the mass and velocity independent parts of the non-zero 
radial and latitudinal components of the inertial forces ( A.8 1— ( A.ll I take the form: 

(a 2 + r 2 )[(r-l)(q 2 + r 2 )-2rA] 4r r 

Gr= ~Kv 7 + (A ' 19) 

Zr = _ 2(fl 2 + r 2 )[(r-l)(a 2 + r 2 )-2rA] _ 4r + r-1 _ 

r Av fi A ' 

2a sin 9 . A . 9 , 9 , nn . 

C r = ^{-rA(3a 2 + 4r 2 + a 2 cos 29) 

[ivy A 

+ 2(a 2 + r 2 )[(2r - l)(a 2 + r 2 ) - a 2 (r - 1) sin 2 6}}, (A.21) 

and 

g g = la 2 S m2e( 4 ---^], (A.22) 
Z e = cot (9 + a 2 sin 2(9 f 2 - — ) , (A.23) 

\fl v ) 

4a 3 cos sin 2 0(A - a 2 - r 2 )\/A 
Ce = ^— , (A.24) 

f i = a 2 + 2r 2 + a 2 cos 29, (A.25) 
i/ = (a 2 + r 2 ) 2 -Aa 2 sin 2 6». (A.26) 

Eulcr force L k vanishes because of the uniformity of the motion. The charge and 
velocity independent parts of the non-zero electric and magnetic force components 



where 
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(A. 14 1 and (A. 15 1 are given by the relations 



£ r = e ( fl2 + r2 )( r2 ~ fl2 cos2 °) ( A 27 ) 



£ 3 / 2 vVA 

aesinOta 2 cos 2 — r 2 ) , , . 

M r = ^ 1, (A.28) 



Be = - 2a2e ;^r 2 ^ (A.29) 



cos 

E 3 /2. 



2aercos6>(Y 2 + a 2 ) , ao „n 
^« = ^TTr: ( A - 3 °) 
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